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(.B* /10 -2 Gauss) 1//2 , where is the value of the primordial field at recom- 
bination. The depolarizing mechanism reduces the damping of anisotropics 
due to photon diffusion on small angular scales. The I ~ 1000 multipoles of 
the CMB temperature anisotropy correlation function in a standard cold dark 
matter cosmology increase by up to 7.5% at frequencies where depolarization 
is significant. 
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1 Introduction 



The cosmic microwave background radiation (CMB) is expected to have acquired a 
small degree of linear polarization through Thomson scattering |]J , which polarizes the 
radiation if there is a quadrupole anisotropy in its distribution function . Typically, 
the CMB degree of linear polarization is expected to be more than ten times smaller 
than the relative temperature anisotropy on comparable angular scales, at least within 
a standard ionization history of the Universe. The CMB has not yet been observed 
to be polarized, the upper limit on its degree of linear polarization on large angular 
scales being P < 6 x 1CT 5 ||. When measured, the CMB polarization will provide 
a wealth of information about the early Universe, additional to that revealed by the 
CMB anisotropy. 

The polarization properties of the CMB may prove particularly valuable to either 
constrain or detect an hypothetical primordial magnetic field @, [|. A cosmological 
magnetic field could leave significant imprints upon the CMB polarization through 
the effect of Faraday rotation. After traversing a distance L in a direction q within 
an homogeneous magnetic field B, linearly polarized radiation has its plane of polar- 
ization rotated an angle 
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n e is the total number- density of electrons and x e its ionized fraction. A is the wave- 
length of the radiation, m is the electron mass, and c is the speed of light. We work 
in Heaviside-Lorentz electromagnetic units (a = e 2 /47r m 1/137 is the fine structure 
constant if we take % = c = 1). 

Faraday rotation of synchrotron emission by distant galaxies serves, for instance, 
to estimate the value of galactic and extragalactic magnetic fields ||. Faraday ro- 
tation acts also as a depolarizing mechanism. If an extended source emits polarized 
radiation, the total outcome may become significantly depolarized by a magnetic 
field, after the radiation emanating from points at different depths within the source 
experience different amounts of Faraday rotation. This process affects significantly 
the radio emission of galaxies and quasars 0. 

In this paper we analyze the depolarizing effect exerted by a primordial magnetic 
field upon the CMB across the last scattering surface. We consider a Robertson- 
Walker universe with scalar, energy-density fluctuations, and assume a standard 
thermal history. We make use of an analytic approach |J, based on a recent re- 
finement and extension |5j of the tight- coupling approximation JlUfl , that highlights 
the physical process responsible for the CMB polarization and its dependence upon 
various cosmological parameters, while still yielding reasonably accurate results. The 
polarization of the CMB is proportional to the width of the last scattering surface 
(LSS), the interval of time during which most of the CMB photons that we observe 
today last-scattered off free electrons. A primordial magnetic field could prevent the 
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polarization from growing across the full width of the LSS. We shall see that the effect 
is controlled by the dimensionless and time-independent parameter 

2vrez/ 2 ' l 10- 3 Gauss^ Vq ) ' ^ 

The coefficient F represents the average Faraday rotation (in radians) between Thom- 
son scatterings Q. z/ is the CMB frequency observed today = B(t*) is the 
strength of the primordial magnetic field at a redshift z* = 1000, around the time 
of decoupling of matter and radiation. Current bounds suggest that a magnetic field 
pervading cosmological distances, if it exists, should have a present strength below 
B « 10~ 9 Gauss ||. It is conceivable that the large scale magnetic fields observed 
in galaxies and clusters have their origin in a primordial field, and several theoretical 
speculations exist about its possible origin [flT]| . A primordial magnetic field is ex- 
pected to scale as B(t) = B (to) a 2 (to)/ a 2 (t), where a(t) is the Robertson- Walker scale 
factor. Thus, a primordial field with strength B* = 10~ 3 Gauss at recombination 
would have a present strength 

A primordial magnetic field may significantly depolarize the CMB right before 
its decoupling from matter. The effect is sensitive to the strength of the magnetic 
field at recombination, not to its present strength. A value of B* somewhat larger 
than 10~ 3 Gauss is not ruled out. Compatibility with big-bang nucleosynthesis, for 
instance, places an upper bound that, extrapolated to the time of recombination, is 



at most B* = 0.1 Gauss ||12jj . Recent proposals for either detecting or constraining a 
primordial field at recombination were suggested in [^|, |T3[ . In |J , Faraday rotation of 
the CMB polarization was analyzed in the limit of small rotation angles, concluding 
that a measurement of the effect could provide evidence for magnetic fields of order 
B* « 10~ 3 Gauss at recombination. In [13|] the change in the photon-baryon sound 
speed in the presence of a magnetic field of order B* = 0.2 Gauss was claimed to 
distort the structure of the acoustic peaks in the CMB anisotropy power spectrum at 
a level detectable by currently planned CMB experiments. 

We shall entertain in our discussions the possibility that the strength of the pri- 
mordial magnetic field at recombination be somewhat larger than B* = 10~ 3 Gauss. 
We will show that currently planned CMB experiments might be sensitive to the effect 
of depolarization upon the temperature anisotropy power spectrum on small angu- 
lar scales if B* is around or larger than 0.01 Gauss, while experiments at somewhat 
lower frequencies would be sensitive to primordial fields of strength around B* « 10~ 3 
Gauss. 

The impact of depolarization upon anisotropy comes about as a consequence of 
the polarization-dependence of Thomson scattering, which feeds back polarization 
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into anisotropy ||, [14], [L5|]. The dominant effect is a reduction in the exponential 
damping due to photon diffusion, which results in an increase of the anisotropy at 
those frequencies for which depolarization is significant. We shall perform an analytic 
estimate of the effect, based on the tight coupling approximation. In order to make 
more quantitative and specific predictions about the impact and potential measura- 
bility of the effect of depolarization upon temperature anisotropy, we shall also use 
a recently developed numerical code [|16|] to integrate the Boltzmann equations in 
a standard cold dark matter model. We shall see that the temperature anisotropy 
correlation function multipoles at / « 1000 increase by up to 7.5% at frequencies 
where depolarization is significant. We conclude that a primordial magnetic field of 
strength around 10~ 2 Gauss at recombination is worth of membership in the list of 
multiple cosmological parameters that one may attempt to determine through CMB 
anisotropy measurements on small angular scales [ |17[| . 

The paper is structured as follows. In section II, we write down and describe the 
radiative transfer equations for the total and polarized photon-distribution function 
in the presence of a single Fourier mode of the scalar metric fluctuations. We include 
the term describing Faraday rotation by a primordial magnetic field. We solve these 
equations in the tight coupling approximation, and find the dependence of the degree 
of polarization upon the frequency of the CMB photons in the presence of a primordial 
magnetic field. In section III we discuss the effects of the depolarizing mechanism 
upon the anisotropy of the CMB on small angular scales, both analytically as well 
as numerically. We discuss the possibility that the effect be detected by currently 
planned CMB experiments. Section IV is the discussion and conclusion. 



2 Depolarization by a magnetic field 
2.1 Boltzmann equations 

We begin by considering the radiative transfer equations for a single Fourier mode 
of the temperature and polarization fluctuations in a Robertson- Walker spatially-flat 
Universe with scalar (energy-density) metric fluctuations, described in terms of the 
gauge-invariant gravitational potentials \1/ and $. We follow the notation and for- 
malism of Ref. H. The total temperature fluctuation is denoted by Ay, while the 
fluctuation in the Stokes parameters Q and U are denoted by Aq and respec- 
tively. The degree of linear polarization is given by Ap = (Aq + A 2 /) 1 / 2 . All three 
quantities are expanded in Legendre polynomials as Ax = + l)Ax ! -P«(A i )5 where 
/i = cos 8 = k ■ q/\k\ is the cosine of the angle between the wave vector of a given 
Fourier mode k, and the direction of photon propagation q. The evolution equations 
for the Fourier mode of wave vector k of the gauge-invariant temperature and polar- 
ization fluctuations || || [18], 19fl , including the Faraday rotation effect of a primordial 
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magnetic field 0, read 

A T + ikfi(A T + = k[A T - A To - fiV b + -P 2 (n)S P ] (4) 

A Q + ik/iA Q = -k[Aq - -(1 - P 2 (/i))5 P ] + 2^ B A l/ (5) 

A[/ + ikjj.Au = -kAu - 2uj b A q . (6) 

We have defined 

S P = -A T2 - A Q2 + Aq (7) 

which acts as the effective source term for the polarization. V b is the bulk velocity of 
the baryons, which verifies the continuity equation 

V b = --V b -iW + ?-(3A TL -V b ) . (8) 
a H 

An overdot means derivative with respect to the conformal time r = / dtao/a, with 
a(t) the scale factor of the spatially flat Robertson- Walker metric, and ao = a (to) its 
value at the present time. R = 3pfc/4p 7 coincides with the scale factor a(t) normalized 
to 3/4 at the time of equal baryon and radiation densities, k = x e n e aTa/ao is the 
Thomson scattering rate, or differential optical depth, with n e the electron number 
density, x e its ionized fraction, and op the Thomson scattering cross-section. Finally, 
u B is the Faraday rotation rate 

dip e 3 n e x e B -da . . 

If there were axial symmetry around k and no Faraday rotation, one could always 
choose a basis for the Stokes parameters such that U — 0. A magnetic field with 
arbitrary orientation breaks the axial symmetry, and Faraday rotation mixes Q and 
U. 



2.2 Tight coupling approximation 

We now solve the equations (0,|,|D in the tight - coupling approximation, which 
amounts to an expansion in powers of krc, where tq = k" 1 is the average conformal 
time between collisions. 

At times earlier than decoupling, Thomson scattering is very efficient, and the 
mean free path of the photons is very short. The lowest order tight-coupling expres- 
sion constitutes in that case an excellent approximation. It implies that the photon 
distribution function is isotropic in the baryon's rest frame, and hence the polariza- 
tion vanishes |§. To first order in krc there is a small quadrupole anisotropy, and 
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thus a small polarization. As decoupling of matter and radiation proceeds, the tight- 
coupling approximation breaks down. Still, for wavelengths longer than the width 
of the last scattering surface, it provides a very accurate approximation to the exact 
result. 

In the absence of a magnetic field (ujb = 0), the tight - coupling solutions, to first 
order in krc, are such that 

A u = ; A Q = ^S P sin 2 9 (10) 

S P = -^A T2 = Ukr c A Tl = ~t c A , (11) 

where we defined A = A To + $. Notice that Aq = — 5Aq 2 = — |A Ta = |Sp, while 
all multipoles with / > 3 vanish to first order in krc- All quantities of interest can be 
expressed, in the tight-coupling approximation, in terms of A , which in turn verifies 
the equation of a forced and damped harmonic oscillator || 

x r R 16 k 2 r c ii k 2 . k 2 ,^ ,„ _ T1 

A » + [lTfl + 45(TTi5) Ao + 5(TT^ A ° = 3(rT^ l *- (1 + fl) * 1 ' (12) 

where we have neglected 0(R 2 ) corrections. 

Now consider the effect of the magnetic field {up ^ 0), assumed spatially homo- 
geneous over the scale of a perturbation with wave-vector k. Faraday rotation breaks 
the axial symmetry around the direction of the wave-vector. The depolarizing effect 
of Faraday rotation depends not only upon the angle between the magnetic field and 
the direction in which the radiation propagates, but also upon the angle between the 
magnetic field and the wavevector k. Nevertheless, we shall only be interested in 
the stochastic superposition of all Fourier modes of the density fluctuations, with a 
Gaussian spectrum that has no privileged direction. Average quantities thus depend 
only upon the angle between the line of sight and the direction of the magnetic field, 
but not upon the angle between the magnetic field and the wavevector k, which is 
integrated away. For simplicity of the calculation, when computing the evolution of 
perturbations with wave-vector k we shall consider a magnetic field with no compo- 
nent perpendicular to k. This choice also satisfies the condition of axial symmetry 
around k, under which eqs. (|,0,||) for At, Aq and Au were derived. We shall later 
use the result of this calculation for the stochastic superposition of all Fourier modes 
with arbitrary orientation relative to the magnetic field. This simplification will result 
at most in an underestimate of the net depolarizing effect, since the case B || k is 
that for which depolarization is less effective, the magnetic field being perpendicular 
to the direction in which polarization is maximum. 

To first order in krc the tight-coupling solutions in the presence of a homogeneous 
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magnetic field B || k are such that: 

4 (1 + F 2 cos 



^ a 3 Sp sin 2 9 

A c/ = -Fcos^A Q ; A Q = -— ^— — — (13) 



where we have defined the coefficient F as 

Fcosfl = 2wbt c (14) 

and so 

e 3 B ( B* w10GHzn2 / 



4ir 2 m 2 (7T v 2 ^10 3 Gauss^ vq 

The coefficient F represents the average Faraday rotation between collisions, since 
2ujb is the Faraday rotation rate and tq = k~ x is the photons mean free path (in con- 
formal time units). When calculating the evolution of each mode k we have assumed 
that the strength of the primordial magnetic field scales as B(t) = B(t*)a 2 (t*)/a 2 (t), 
which is justified by flux conservation and because the Universe behaves as a good 
conductor 0. Since the frequency also redshifts as v = v a(to) / a(t) , the parameter F 
is time- independent. Uq is the frequency of the CMB photons at present time, while 
B* is the strength of the magnetic field at a redshift = 1000, around recombina- 
tion. Within a standard thermal history, with no early reionization, depolarization 
is only significant across the LSS, and it thus depends only upon the value of the 
primordial magnetic field around the time of recombination. Notice that Faraday 
rotation between collisions becomes considerably large, paving the way to an efficient 
depolarizing mechanism, at frequencies around and below v& defined such that 



so that 



u d « 8.4 GHz ( ^ ) 1/2 ^27GHz( ^ Y /2 . (17) 

Vio-3Gauss^ Vo.01 Gauss' v ; 

From eqs. (|i~3|) we can read the values of Aq and Aq 2 . They reduce to eqs. (|T0D 
with 0(F 2 ) corrections for small F, while they vanish as F~ l for large F. We write 
them as: 

A Qo = \do{F)S P ■ A Q2 = -±d 2 (F)S P . (18) 

The coefficients do,d 2 are defined so that di ~ 1 + 0(F 2 ) for small F, while di 
0(1/ F) as F — > oo, and represent the effect of depolarization. They read: 

w = 8 [ «^ot (1 + ^ ) _^ (19) 

dj(F) = ^[fE^CQ (1 + ^ + ^ ) _^_^ . (20) 
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In terms of the combination 



5., d 2 > 15 rarctan(F) 2 1 5 1 n 

and using the definition of Sp, we find the relation 

A T2 = -S P (l - |d) (22) 
5 

and from the equation for Ap in the tight coupling limit we get 

4 4 

Sp = —, -ikTcATi = ; —TcAn . (23) 

3(3 -2d) 1 3(3 -2d) 60 v ; 

Notice that d « 1 — F 2 /7 if F << 1 while d — > jf^F -1 for large F. We stress 
here again that in a general case the depolarizing coefficient d depends upon the 
angle between k and B. The net anisotropy and polarization being the outcome of 
the stochastic superposition of all Fourier modes of the density-fluctuations, with 
a spectrum that has no privileged direction, the average depolarizing factor, after 
superposition of all wavevectors in arbitrary orientations with respect to the magnetic 
field, depends only upon F. The average depolarizing factor might slightly differ from 
that calculated with k || B, which at most underestimates the average effect. 



Equations ( |T3|j22]j23| ) condense the main effects of a magnetic field upon polariza- 
tion. When there is no magnetic field (F — 0, d — 1) Au = and Aq = — yAj- 2 sin 2 6. 
A magnetic field generates Ay, through Faraday rotation, and reduces Aq. In the 
limit of very large F (large Faraday rotation between collisions) the polarization van- 
ishes. The quadrupole anisotropy A Tz is also reduced by the depolarizing effect of 
the magnetic field, by a factor 5/6 in the large F limit, because of the feedback of 
Aq upon the anisotropy or, in other words, because of the polarization dependence 
of Thomson scattering. The dipole A^ and monopole Ar are affected by the mag- 
netic field only through its incidence upon the damping mechanism due to photon 
diffusion for small wavelengths, that we shall discuss in detail in section III. Indeed, 
the equation for A = A To + neglecting 0(R 2 ) contributions, now reads: 

x r R 16 (5 -3d) k 2 r c 1A k 2 . k 2 ^ ,„ , ^ 

(24) 

The damping term is reduced by a factor 5/6 at frequencies such that d << 1, for 
which depolarization is significant. 

We have assumed that the magnetic field is spatially homogeneous. We can expect 
corrections to our result if the field is inhomogeneous over scales smaller than tq 
at any given time around decoupling. Indeed, if the field reverses its direction N 
times along a photon path during a time tq, Faraday rotation will not accumulate 
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as assumed above. In that case depolarization would start to be significant only at 
those frequencies such that Faraday rotation is large over the scale on which the 
magnetic field reverses its direction. The frequencies at which depolarization starts 
to be significant would thus be reduced by a factor l/yN. 



2.3 Frequency-dependence of the degree of polarization 

The anisotropy and polarization observed at present time can be evaluated using the 
formal solutions of eqs. (f|,§,§]) 



(25) 



At(t ) = f T0 dre ik ^-^g(T)[A To (r) + »V b (r) - |P 2 (/i)S P (r)] 
+ J TO dTe ikfl( - T - To) e- K( - TO ' T) (4> - $) 

A Q (r ) = jf dre^-^g{r){^[l - P 2 (fi)]S P (r) + FA^(r)} (26) 

Au(t ) = - r dre lk ^-^g(r)FA Q (r) (27) 
Jo 

where 

g(r) = ke- K{T(hT) (28) 

is the visibility function. It represents the probability that a photon observed at r 
last-scattered within dr of a given r. For a standard thermal history, with no signifi- 
cant early reionization after recombination, g(z) is well approximated by a Gaussian 
centered at a redshift of about z ~ 1000 and width Az ~ 80 p0| . In conformal time, 



we shall denote the center and width of the Gaussian which approximately describes 
the process of decoupling by To and Atd respectively. 

The visibility function being strongly peaked around the time of decoupling, the 
first integral in eq. fl25|) for the anisotropy is well approximated, at least for wave- 
lengths longer than the width of the last scattering surface, by its instantaneous 
recombination limit. In that case it reduces to the tight-coupling expression of its 
integrand evaluated at time r = m @- This first integral is dominated by its first 
two terms, proportional to the monopole At and the baryon velocity Vb (in turn 
proportional to A Pl ) respectively. The quadrupole term Sp gives a negligible contri- 
bution for long wavelengths, but becomes relatively significant on small scales. The 
second integral in eq. ( j25"D corresponds to the anisotropies induced by time-dependent 
potentials after the time of last-scattering. 

Eqs. (^,^) for the polarization can be approximated replacing the integrand by 
its tight-coupling expression. Then 

A O (r ) = 2m / dre^ e ^g(r)S P (r) (29) 

4 (1 + F l cos" 1 0) Jo 



8 



A u (r ) = -Fcos6A Q (r ) (30) 
while the total polarization, Ap = (Aq + A 2 /) 1 / 2 reads 

A P (r ) = Vl + F 2 cos 2 flA Q (r ) (31) 

Evaluation of the time integral in eq. (|29|) requires a more detailed knowledge of 
the time-dependence of the integrand than in the case of the anisotropy. Indeed, 
the tight-coupling expression (p3|) for the quadrupole term Sp being proportional to 
the mean free path tc, which varies rapidly during decoupling, the instantaneous 
recombination approximation becomes inappropriate. The induced polarization is, 
indeed, proportional to the width of the last scattering surface. Adapting the method 
of |J to include also the effect of the primordial magnetic field, we write down the 
equation satisfied by Sp when all other quantities are already approximated by their 
first-order tight-coupling expressions: 

S P + — (3 - 2d)kS P = -ikA Tl (32) 
10 5 



S P (r) = \ik f ' dr'A Tl e-^M(^ (33) 



Neglect of Sp returns the tight-coupling result of eq. (p3|). Instead, the formal solution 
to equation ( |32"1) 

> rr' 

Ik I dT* At. P. 

5 

tracks down the time-dependence of Sp through the decoupling process with better 
accuracy. 

For wavelengths longer than the width of the LSS we can neglect the time variation 
of Ap x and that of e lfccos6, ( r - r o) around decoupling. We also approximate the visibility 
function by a Gaussian, which justifies the approximation k(tq,t) ~ ~ K ^^ ||21|| . 
Then 

Sp{T)^-ikA Tl {T D )Ar D e^^- 2 ^ [°° ^ e -^~^) (34) 
5 Ji x 

k(tq,t) 

approximations, 



where the integration variable has been changed to x — K y°' T> , Thus, within these 



£ drg(T)Sp(r) = -fifcA Tl (t d )At d J q °° ^e"^ K Ei(- A(3 - 2d) K ) 

(35) 

= ^A Tl (r D )Ar D [ln(f ) - ln(3 - 2d)} 

Finally, the total polarization induced at an angle 9 with respect to the wavector 
k, reads 

o in n ik cos e(m -to) 

A ^-(TT65)H T )-M3-2 rf) ] V1 + FW , M T ,MA, . (36) 
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It can also be written as follows, in terms of the polarization that would be induced 
if there were no magnetic field (or equivalently, in terms of the polarization at fre- 
quencies large enough such that the depolarizing effect is negligible): 

A P (0, F) = D(6, F)A P (B = 0) (37) 

where we have defined the depolarizing factor as 

D(6, F) = 1 =f(F) (38) 

V ; VI + F* cosV V ; V ; 

with 

f(F) - 7 h M3-2d) 1 

n ^~l + 6d^ ln(10/3) J [ } 

Eq. Q3"8D, together with the defining eqs. ( |T5| , |r6[ ) and ( ]2~I| ) for F and d, summarize 



the main result of this section. Notice that 

/ -> 1 as F -> ; / -> 0.61 as F -> oo (40) 

The polarization observed at present times depends upon the angle between the 
line of sight and the orientation of the magnetic field at the time of decoupling. There 
will be no depolarization if the magnetic field is perpendicular to the line of sight. 
The magnetic field is likely to change orientation randomly over scales longer than 
the Hubble radius at the time of decoupling, which subtends an angle of order one 
degree in the sky, so that after averaging over many regions separated by more than 
a few degrees, we can always expect a net average depolarizing effect. To roughly 
estimate its order of magnitude we could assume an average component of B parallel 
to the observation direction of order B / v2 and define an average D as 

D= i 1 = f(F) . (41) 

Fig. 1 displays the depolarizing factor D as a function of the CMB frequency z/ . We 
have plotted it for three different values of the magnetic field B* to help visualize 
the relevant frequency range, but notice that since depolarization depends only upon 
F = (z/d/z/ ) 2 , the plot for an arbitrary value of B* is identical to that corresponding 
to another value of the magnetic field after an appropriate scaling of the frequency 
units, proportional to the square root of the magnetic field. 

At low frequencies, those for which the effect is large, the average depolarizing 
factor scales as 

D w 0.6— w 0.85f— ) 2 if u « u d . (42) 
F v zV 

At comparatively large frequencies instead 

D fx 1 — 0.36-F 2 = 1 — 0.36(— ) 4 if u » v A . (43) 
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3 Effects upon the anisotropy 



Depolarization by a primordial magnetic field has significant and potentially measur- 
able effects upon the anisotropy of the CMB on small angular scales. Indeed, the 
polarization properties of the CMB feed back into its anisotropy, as evidenced in eq. 
(^), due to the polarization dependence of Thomson scattering. The dominant effect 
of polarization upon anisotropy derives from its impact upon the photon diffusion 
length H [TJ], |0J, which damps anisotropies on small angular scales || E3I. It was 
shown in JL4] , through numerical integration of the Boltzmann equations, that neglect 
of the polarization properties of the CMB leads to an overestimate of its anisotropy 
on small angular scales as large as 10%. We thus expect depolarization by a primor- 
dial magnetic field to introduce a significant frequency- dependent distortion of the 
CMB anisotropy power spectrum. Notice that a different (frequency-independent) 
distortion of the CMB anisotropy power spectrum by a primordial magnetic field, 
due to its impact upon the photon-baryon fluid sound speed, was recently discussed 
in HI. 



3.1 Reduced diffusion-damping 

Photon diffusion damps anisotropies on small angular scales ||, |22|, ^3j. The effect 
is described by the term proportional to k 2 rc Ao in eq. (|24]) . The photon diffusion 
length depends upon the degree of polarization of the CMB @, [14], Thus, the 
photon-diffusion length is different at frequencies where the depolarizing effect is 
significant. 

The damping of anisotropies on small angular scales due to photon diffusion can be 
found, now including the full independence, by solving the tight-coupling equations 
to second order, assuming solutions of the form 

A x (r) = A x e^ T (44) 

for X = T, Q, and U, and similarly for the baryon velocity V&. One then finds that 

k 

u = + z 7 (45) 

'3(1 + R) 



with the photon-diffusion damping length-scale determined by 

k2 k2r c ( § (5 -3d) R 2 s 



The depolarizing effect of the magnetic field reduces the viscous damping of anisotro- 
pies. In the case of small R, such that R 2 terms can be neglected, the damping factor 
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7 is smaller by a factor 5/6 at those frequencies for which the depolarizing effect is 
large. 

We make now an analytic estimate of the effect of the frequency-dependence of 
the photon diffusion length, in the presence of a primordial magnetic field, upon the 
CMB anisotropy power spectrum. The temperature anisotropy correlation function 
is typically expanded in Legendre polynomials as 

1 oo 

C{6) =< A T (ni)A T (n 2 ) >n 1 -n 2 = cos e= -r- £(2Z + l)Ci^(cos0) • (47) 

47r 1=0 

The multipole coefficients of the anisotropy power spectrum are given by 

Q = (4tt) 2 | k 2 dkP(k)\A Tl (k,T Q )\ 2 (48) 

with P(k) the power-spectrum of the scalar fluctuations, assumed scale-invariant in 
the sCDM model. The largest contribution to a given multipole C\ comes from those 
wavelenghts such that / = k(r — To), where r is the conformal time at present and rp, 
the conformal time at decoupling. The average damping factor due to photon diffusion 
upon the C/'s is given by an integral of e -27 times the visibility function across the 
last scattering surface || . It depends upon cosmological parameters, notably R, and 
upon the recombination history. Approximately, and for a standard cold dark matter 
model, we can take 2 / y(d = 1) « (//1500) 2 . The relative change in the C;'s due to 
the change in the photon-diffusion length, as we move down from frequencies where 
depolarization is insignificant (d = 1) to lower frequencies (d « 1), is then given by 

In figure 2 we plot ACi (expressed as a percentage) at I = 1000 as a function of 
frequency, for three different values of the magnetic field = 0.001, 0.01 and 0.1 
Gauss. Once again, the graph for an arbitrary value of 5* can be read from any of 
these with an appropriate scaling of the frequency units. We have chosen to display 
the effect at I = 1000, that will be accessible by the recently funded CMB satellite 
experiments, MAP [51 and COBRAS /SAMBA 



3.2 Reduced quadrupole contribution 

The depolarizing effect of a primordial magnetic field also changes the strength of the 
quadrupole term Sp around decoupling, and thus its incidence upon the anisotropy 
of the CMB on small angular scales. Indeed, the quadrupole anisotropy and the 
polarization of the CMB at the time of recombination contribute to the presently 
observed anisotropy through the following term of eq. ( f25j ) 

A Sp (r ) = -ip 2 (cos#) r dre ikcos0 ^g(r)Sp(r) . (50) 
Z Jo 
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This term is negligible for long wavelengths, those that dominate the lowest multipoles 
of the present anisotropy, but becomes non-negligible on small angular scales (large 
multipoles). Indeed, in the tight coupling approximation Sp oc TcA and thus, barring 
a very strong time-dependence of the scalar potential, the contribution of Sp is well 
below that of the monopole term, except for small wavelengths. 

The depolarizing effect of a magnetic field modifies the value of Sp around de- 
coupling, compared to what it would have been if there were no magnetic field, and 
then 

A Sp (F) = f(F)A Sp (B = 0) (51) 

with / as defined in eq. QBTf). Recall that / ~ 1 if v » i/ d while / ~ 0.6 if 
v << v&. Thus, at frequencies such that the depolarizing effect of the magnetic field is 
significant, the partial contribution of the quadrupole term Sp to the total anisotropy 
is reduced by a factor 0.6 compared to that at frequencies where depolarization is 
unimportant. On small angular scales this could represent a decrease of the anisotropy 
by a few per cent. The effect is opposite to that of the change in diffusion damping, 
but is likely to be less significant on small angular scales. 



3.3 Numerical estimate of the effect upon the anisotropy 

In order to accurately ascertain the net effect of the depolarizing mechanism upon 
the CMB anisotropy and to make definite quantitative predictions within a standard 
cosmological model, we turn now to the numerical integration of the Boltzmann equa- 
tions (|4]|^,|). We use the recently developed code CMBFASTflBH, that integrates the 



sources over the photon past light cone. Its starting point are the formal solutions 



( |25| , P6|j27D , where the geometrical and dynamical contributions are separately handled 
to improve efficiency As in our analytic estimates, when computing the evolution of 
each Fourier mode we introduce in the code the Faraday rotation term with the angu- 
lar dependence corresponding to the case where B has no component perpendicular 
to k. 

Figures 3 and 4 summarize the numerical calculation of the effect of depolarization 
upon temperature anisotropy, in a standard cold dark matter model (sCDM). 

The quantity plotted in Fig. 3 is 1(1 + l)Cj, for the sCDM model without a 
magnetic field and with a magnetic field and at frequencies such that F = 1,4,9, 
corresponding to u = v&,v&/2 and i>d/3 respectively. Fig. 3 clearly shows that the 
CMB anisotropy on small angular scales increases at frequencies where depolarization 
is significant. This result indicates that the reduction in diffusion damping due to 
depolarization is the dominant effect among the two opposite effects discussed in the 
previous sections. 

Fig. 4 displays the same results but expressed in terms of AC;, the percentual 
increase in C\ relative to the case without magnetic field. The monotonic curves 
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in the same figure, included for comparison purposes, correspond to the analytic 
estimate of the effect of reduced diffusion damping, eq.(f|^). As expected, the effect is 
larger on smaller angular scales (larger I) . The numerical result approximately follows 
the analytic estimate of the effect of the reduction in diffusion damping. The total 
effect, however, does not increase monotonically with I. This can be understood as 
a consequence of the nature of the subdominant quadrupole contribution Sp, which 
oscillates out of phase with the Cz's || (remember that Sp oc A ), and is reduced by 
depolarization through a factor f(F). 

It is also clear from Fig. 4 that the analytic result for the change in diffusion 
damping due to depolarization overestimates the total effect at high 1. This is because 
the actual damping in the G\ spectra has two contributions, one from Silk damping 
and the other due to cancellations in the integral across the last scattering surface 
produced by the oscillations in the exponential and sources in equation ( p5f) . Only 
Silk damping is reduced by the magnetic field, and that is why equation (|49|) slightly 
overestimates the net effect. 

The analytic and numeric calculations are in very good agreement around I ~ 
1000. The frequency-dependence of ACi at I = 1000 plotted in Fig. 2 fits very well 
the analogous result after the full numerical integration of the Boltzmann equations. 

We conclude that the depolarizing effect of the magnetic field results in an increase 
of the anisotropy correlation function multipoles of up to 7.5% (for sufficiently low 
frequencies) on small angular scales (I ~ 1000), those that will be accessible by future 
CMB satellite experiments such as MAP and COBRAS/SAMBA. The frequencies at 
which the effect is significant, however, depend on the strength and coherence length 
of the primordial magnetic field at the time of recombination. 

Depolarization depending upon frequency, the effect might be difficult to separate 
from foreground contamination. The relative change of the Cj's at / = 1000 is larger 
than 2% on frequencies below 30 GHz (accessible to the first two channels in MAP), 
if .B* = 0.02 Gauss or larger. The first two channels in COBRAS/SAMBA being 
at 31.5 and 53 GHz, the signal would reach a 2% level within this range only if B* 
is around or larger than 0.1 Gauss. However, COBRAS/SAMBA might reach out 
to larger values of I, where AC/ is larger, and might thus have a sensitivity to the 
depolarizing effect of B* comparable to MAP. In any case, both experiments will be 
sensitive to a magnetic field around B* =0.1 Gauss, and would thus at least be able 
to place a direct constraint on B* comparable or better than the one obtained from 
extrapolation of the nucleosynthesis bound. 

Experiments searching CMB anisotropy and polarization at smaller frequencies, 



which currently operate down to 5 GHz [26, 071, may play a significant role to detect 



the depolarizing effect of a primordial magnetic field. 
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4 Conclusion 



The CMB is expected to have a small degree of linear polarization. Several esti- 
mates were made for the predicted polarization, both in the context of anisotropic 
cosmological models [[l], j28fl , as well as in isotropic and homogeneous cosmologies per- 
turbed with either energy- density fluctuations or gravitational waves 2T , |2TS| | . The 



1-5 



polarization of the CMB remains undetected, its upper limit being P < 6 x 10 

A primordial magnetic field depolarizes the CMB radiation on those frequencies 
that experience a significant amount of Faraday rotation around the time of decou- 
pling. In this paper, we have applied the analytic method developed in Ref. || 
to estimate the depolarizing effect of a primordial magnetic field across the last- 
scattering surface, assuming a standard ionization history. The result is expressed by 
eqs. ([!^,[4^,[43|) and is represented in Fig. 1. The CMB becomes significantly depo- 
larized at frequencies around and below 30 GHz (B*/0.01 Gauss) 1//2 , below which the 
degree of polarization decreases quadratically with frequency. B* is the value of the 
primordial field at a redshift = 1000, around recombination, likely to be 10 6 times 
larger than an hypothetical cosmological magnetic field at present times. 

The average depolarizing factor depends only upon the parameter F, as defined by 
eq. (15), which represents the average Faraday rotation between collisions. We have 



calculated the depolarizing factor d(F), as given by equation fl2~T|), in the particular 
case of a wavevector k || B. In a general case, the factor d depends upon the angle 
between k and B. This dependence integrates away in average quantities, after the 
stochastic superposition of all Fourier modes of the density fluctuations. The value 
derived here for d is at most an underestimate of the average depolarizing effect, 
which would eventually start to be significant at slightly larger frequencies. Our 
derivation also assumed that Faraday rotation accumulates over the width of the last 
scattering surface. If the primordial magnetic field is very entangled over that scale, 
the depolarizing effect starts to be significant at smaller frequencies. 

The depolarizing mechanism has a significant effect upon the anisotropy of the 
CMB on small angular scales. On those angular scales and at frequencies such that 
the depolarizing effect is large, the damping of anisotropies by photon diffusion is 
reduced, which results in a significant increase of the anisotropy at a fixed angular 
scale. Besides, depolarization reduces the contribution of the intrinsic quadrupole 
anisotropy. Figure 2 displays the estimate for the percentual change of the anisotropy 
power spectrum at / = 1000 due to the reduction in diffusion damping, as a function 
of frequency and for different values of the primordial magnetic field at recombination. 

We conclude that a primordial magnetic field increases the anisotropy of the CMB 
by up to 7.5% at I ~ 1000 in a standard CDM cosmology. The asymptotic strength of 
the effect is independent of the intensity of the magnetic field, but the frequencies at 
which it starts to be significant are those around and below 30 GHz (.B^/O.OlGauss) 1 / 2 . 
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Measurements of anisotropy and polarization at sufficiently low frequencies could 
probe primordial magnetic fields in an interesting range. 
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Figure 1: The average depolarizing factor D as a function of the CMB frequency vq. 
The corresponding figure for an arbitrary value of B* is identical to any of these after 
a scaling of the frequency units, proportional to B^ 2 . 
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Figure 2: Analytic estimate of the percentual change due to reduction in diffusion 
damping of the / = fOOO anisotropy correlation function multipoles as a function 
of the CMB frequency for different strengths of the primordial magnetic field at 
recombination. The corresponding figure for arbitrary can be obtained from any 
of these after a scaling of the frequency units, proportional to B X J 2 
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Figure 3: Numerical integration for the multipoles of the anisotropy correlation func- 
tion in a standard CDM model without a primordial magnetic field (F = 0), and 
with F — 1, 4, 9, which correspond to u = u d , u A /2, i/ d /3 respectively, with u A w 
27 GHz (B./O.OlGauss) 1 / 2 . 



21 




Figure 4: Numerical result for the percentual change of C\ as a function of / relative 
to its value without magnetic field in a standard CDM model. The monotonic curves 
also shown for comparison purposes correspond to the analytic estimate of the effect 
of reduced diffusion damping. 
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